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ABSTRACT 

HD 98088 is a synchronised, double-lined spectroscopic binary system with a mag- 
netic Ap primary component and an Am secondary component. We study this rare 
system using high-resolution MuSiCoS spectropolarimetric data, to gain insight into 
the effect of binarity on the origin of stellar magnetism and the formation of chemical 
peculiarities in A-type stars. Using a new collection of 29 high-resolution Stokes VQU 
spectra we re-derive the orbital and stellar physical parameters and conduct the first 
disentangling of spectroscopic observations of the system to conduct spectral analysis 
of the individual stellar components. From this analysis we determine the projected 
rotational velocities of the stars and conduct a detailed chemical abundance analysis of 
each component using both the Synth3 and Zeeman spectrum synthesis codes. The 
surface abundances of the primary component are typical of a cool Ap star, while those 
of the secondary component are typical of an Am star. We present the first magnetic 
analysis of both components using modern data. Using Least-Squares Deconvolution, 
we extract the longitudinal magnetic field strength of the primary component, which is 
observed to vary between -1-1170 and -920 G with a period consistent with the orbital 
period. There is no field detected in the secondary component. The magnetic field in 
the primary is predominantly dipolar, with the positive pole oriented approximately 
towards the secondary. 

Key words: stars: magnetic fields, stars: abundances, stars: chemically peculiar, 
(stars:) binaries: spectroscopic, stars: individual: HD 98088 



1 INTRODUCTION 

Ap stars are magnetic intermediate-mass main sequence 
stars that exhibit distinctive chemical peculiarities in their 
atmospheres. Although other classes of chemically peculiar 
stars exist (e.g. Am stars, Hg-Mn stars), these stars have 
been demonstrated to lack strong, organi s ed magnetic fields 
at their surfaces (e.g. fshorlin et al.ll2002l : IWade et al.]|2006l : 



iMakaganiuk et all 12011 ). Ap stars appear to be the only 
class of middle main-sequence stars fo r which, in all cases , 
an observable magnetic field is present (|Auriere et al.|[2"007t ) . 
The magnetic fields of Ap stars have been established 
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to have important global dipole components ranging in 
strength from hundreds to tens of thousands of gauss. The 
symmetry axis of the dipole component is almost always 
tilted relative the stellar rotation axis. In addition, Ap stars 
generally spin much more slowly than non-peculiar stars of 
similar masses (Abt & Morrcll 1995; Stcpich 2000), and as 
they spin they exhibit line profile variations attributed to 
rotational modulation of patchy, non-axisymmetric lateral 
and vertical distributions of chemical abundance in their 
photospheres. While it is generally accepted that the funda- 
mental mechanism responsible for the chemical pe culiarities 
is mi croscopic chemical diffusion (as described by iMichaudI 
Il970l ). the origin of chemical patchiness, and the relationship 
to the magnetic field, is poorly understood. 



Although mo dern surveys l|Gerbaldi et al.l Il985l : 
[Carrier et al.l [2002') indicate that Ap stars exhibit about 
the same binary frequen cy as normal A-type stars (47 ±5%, 
iJaschek fc G6mej[l970l ). they generally agree that there 
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is a clear lack of Ap sta rs in short-pe r iod (i .e. Porb = 3.0 
days or less) systems (| Carrier et al.l l2002h . One might 
hypothesize that this could be because these systems tend 
to be synchronized. Such a short-period binary would 
imply a short rotational period for both components. Rapid 
rotation would induce rotational mixing, inhibiting the 
development of chemical peculiarities. However, no similar 
absence of short period systems exists for Am binaries. This 
suggests that binaries with short periods are capable of 
hosting chemically peculiar stars, but contain few magnetic 
(chemically peculiar) stars. The few magnetic Ap stars 
reported to be in short period binaries could therefore 
be very interesting targets for studying the origin and 
evolution of magnetic fields in intermediate-mass stars. 

Systems containing two A-type stars are of particu- 
lar interest, since they provide an opportunity for study- 
ing not only the incidence and evolution of stellar mag- 
netic fields, but also of chemical peculiarity in both stellar 
components. A comprehensive literature review reveals that 
these ideal systems are i ndeed very rare. Both HP 135728 
iFrevhammer et al.ll2008h and HD 59435 (|Wade et al.lll999l ') 
have magnetic Ap secondary components; however their pri- 
mary components are giants of spectral type G8. While these 
primary giant G8 stars evolved from main sequence stars 
with masses similar to that of the Ap star, the evolved pri- 
mary component makes these systems less useful to study. 
There are only five reported SB2 binaries with a magnetic 
Ap primary component and a main sequence A star com- 
panion. T he first of these systems is HD 55719, which was 
studied bv lBonsackll|l976h . although as a southern object it 
has received relatively less attention in the modern litera- 
ture HD 5550, HD 22128, and HD 56495 (e.g. lCarrier et al.1 
|2002| ) show composite A-type spectra, although the mag- 
netic nature of the Ap component has yet to be definitively 
established. The final system is HD 98088, which is the sub- 
ject of this paper. This system is very well-suited for the 
study of magnetic fields and chemical peculiarities, and is 
one of the brightest and best-studied magnetic binary sys- 
tems. 

HD 98088 (HR 4369, HIP 55106 ) was first observed as 
a spectroscopic binary bv lAbtl (|l953l '). although he reported 
it to be an SBl. The primary radial velocity was found to 
vary with an orbital period of 5.905 ± 0.005 days. Spec- 
trophotometric analysis using the Mills spectrograph showed 
variability in b eryllium and titanium lines. Five years later 
lBabcocl3 (|l958f ) confirmed the presence of a magnetic field 
in HD 98088 using polarized spectroscopy. The strength of 
the longitudinal magnetic field was found to vary between 
— 1000 and +800 gauss with a peri od consistent with the 
orbital period derived bv I Abj l| 19531 ). These measurements 
made HD 98088 the first observed synch ronized spectro - 
scopic bina ry containing a magnetic star (IBabcoc H i 19581 ). 



In his 1 19651 paper on spectral classification. lOsawal classified 
the Ap primary star to be of spectral subtype SrCr. 

The first complet e orbital solution of this SB2 sys - 
tem was computed bv lAbt et all l|l968t ). lAbt et al.1 (|l968l ') 
searched for photometric eclipses without success. From the 
the difference in equivalent widths of the Ca ll K, Mg ll 
A4481, Na l Di and D2, Si 11 A6371, and Ha lines, the mag- 
nitude difference of the components were estimated to be 
Amv = 1.2. This led them to conclude that the primary 
has spectral type A3Vp and the secondary has spectral type 



A8V. The masses of each component were then inferred from 
the spectral type, and combined with the dynamic Msiv?i to 
determine the system's incl ination angle i = 67°. An inter- 
esting feature discovered bv lAbt et al.1 (| l968l) is that despite 
its short orbital period, HD 98088 has significant orbital 
eccentricity. This suggests that apsidal motion should take 
place. However. I WolfB (i 19741 ) obtained new spectroscopic ob- 
servations that were used in conjunction with older measure- 
ments to refine the period of apsidal motion to at least 700 
years. 



[Carrier et al. 

1 (|200^ reobserved the system using the 
CORAVEL velocimeter. With these measurements, and all 
other published radial velocity measurements of the system, 
the orbital period of HD 980 88 was refined to 5.905111 ± 
0.000004 days. The [Carrier et al. (2002) solut ion is consis- 
tent with the solution of Abt et al.l l|l968i ). [Carrier et al.l 
( 200^ ) observed no significant variation of the longitude of 
periastron, and concluded that the apsidal period is longer 
than one millennium. They used the Hipparcos distance and 
effective temperatures determined from Geneva photome- 
try to derived the masses of the components from their 
HR diagram positions to be Mi = 2.261 ± 0.093 Mq and 
M2 = 1.677±0.085 Mq. Assuming negligible Zeeman broad- 
ening and a synchronized system, they found the radii to be 
Ri = 3.27 ± 0.43 Rq and R2 = 2.10 ± 0.28 Rq using the 
projected rotational velocities and the orbital period. The 
inclination angle was found to be i = 66° by comparing the 
minimum dynamical masses to the evolutionary masses. 

The HD 98088 system potentially provides a very in- 
teresting and rare insight into chemical peculiarities and 
magnetism in A type stars, in the context of close binary 
systems. The similarities between the physical parameters 
of the primary and secondary make it an excellent labo- 
ratory for comparative astrophysics. In this paper we use 
spectroscopic and polarimetric data to improve the orbital 
parameters of the system, to perform a detailed abundances 
analysis of the spectra of both stellar components, and to 
analyze the magnetic properties of the system. 



2 OBSERVATIONS 

Observations of HD 98088 were obtained with the Multi- 
Site Continuous Spectroscopy (MuSiCoS) spectropolarime- 
ter, a now-decommissioned high resolution spectropolarime- 
ter located at the the Telescope Bernard Lyot at the Ob- 
servatoire du Pic du Midi, France. The instrument consists 
of a bench mounted cros s-dispersed echelle spectrograph 
IIBaudrand fc Bohml I1992I ) fibre-fed from a Cassegrain 
mounted polarimeter unit l|Donati et al.|[l999l ). The instru- 
ment provides continuous wavelength coverage from about 
4500 to 6500 A at a resolution of = 35 000. Observations 
were obtained in spectropolarimetric mode, providing Stokes 
V , Q and U spectra as well as S tokes / spectra. T he data 
were reduced using the ESpRIT (iDonati et al.lll997l ) reduc- 
tion tool, which performs relevant calibrations and optimal 
ID spectrum extraction. 

A total of 29 spectra - 10 Stokes V, 10 Stokes Q and 9 
Stokes U - were obtained between Jan 1999 and Jan 2002. 
A summary of the observations is presented in Table [T] 
The typical signal-to-noise ratio (S/N) per 4.5 km/s spec- 
tral pixel was about 170 for Stokes V and 230 for Stokes 
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Q/U . The orbital phases are calculated with the ephemeris 
described in Tabled 

Least-Squares Deconvolution (LSD) was applied to each 
of the calibrated and normalised spectra. Line masks were 
constructed using the derived atmospheric parameters (de- 
termined in Sect.fS} and abundances (derived in Sect.[S]) of 
both the primary and secondary components. Due to the 
similarity between the temperatures and abundances of the 
two stars, only small differences were observed in the profiles 
obtained from the two lines masks. Ultimately, we decided 
to use only the line mask corresponding to the primary to 
extract the LSD profiles. 

All LSD profiles are normalised to a wavelength of 525 
nm and a lande factor of f .41, corresponding approximately 
to the S/N- weighted mean of the observed lines included in 
the mask. 



Table 2. Best fit orbit al parameters, allow ing for offsets between 
the radial velocities of ICarrier et al.l 1I2OO2 ) and older measure- 
ments (Ai)i), and between our measurements and older values 
(Al,2). 



Reduced 


1.2039 


RMS (A) 


2.1863 


RMS(B) 


6.9670 


P(d) 


5.9051102 ± 0.0000023 


To (HJD -2,400,000) 


34401.369 ± 0.015 


Ka (kms-i) 


73.31 ±0.26 


Kb (kms-i) 


100.70 ±0.43 


Vb (kms-i) 


-8.65 ±0.29 


e 


0.1840 ±0.0025 




313.25 ±0.89 


Avi (kms ^) 


-1.34 ±0.38 


Av2 (kms~^) 


-3.02 ±0.40 


Ma /Mb 


1.374 ±0.008 


MAsin^i (Mq) 


1.772 ±0.013 


Mb sin^ i (Mq) 


1.2897 ±0.0089 


a\ s'mi (Rq) 


8.411 ±0.030 


Ob sini (Rq) 


11.553 ±0.050 



3 ORBIT 



4 SPECTRAL DISENTANGLING 



Orbital parameters for HD 98088 were determined using new 
and previously published measurements of the radial veloc- 
ities of both the prim ary and sec o ndary com p onent s. Pub- 
lished data are from IaM (Il953l ). lAbt et all (|l968l ). IWolfj 
l| 19741 ) ■ and [Carrier et al.1 (|2002l ). The new radial velocities 
measurements were obtained using the LSD profiles dis- 
cussed in the previous section, extracted from each of the 
reduced spectra. To measure the radial velocities of the LSD 
line profile s of each componen t we used an IDL code, de- 
scribed bv lAlecian et al.l (|2008l ). that fits the binary profile 
by least-squares using the IDL function CURVEFIT. The 
binary profile is assumed to be the sum of two photospheric 
functions, one for each component of the system. Each pho- 
tospheric function is the convolution of a rotation profile 
and a Gauss ian of instrumental and turbulent velocity width 
ICr^liooi). For the purpose of the fit, the turbulent veloc- 
ity was fixed at 2 kms~^, and the rotation profile broaden- 
ing was left as a free parameter. The errors were calculated 
using a Monte-Carlo-type approach in which five thousand 
observed profiles were calculated by adding random noise to 
the best fit profile. The fitting procedure was repeated on 
each of these five thousand profiles and the standard devi- 
ation of these measurements was adopted as the error. The 
derived radial velocities are reported in Table [T] 

The full set of radial velocity measurements was used 
to determine orbital parameters for the binary system. This 
was done by fitting orbital radial velocity curves to the mea- 
sure d velocities. A small constant offset between our values 
and [Carrier et al.l (|2002l ) was used (~ 2 kms~^), signifi- 
cantly impr oving our fit to the radial velocities, following 
the work of [Carrier et al. who found it necessary to 

shift their radial velocities with respect to previous measure- 
ments. The best fit orbital parameters are presented in Table 
[2I and the fit to the radial velocities is shown in Fig. [T] Our 
orbita l parameters are consistent with those of lCarrier et al.l 
l|2002l ). although our parameters are more precise due to our 
larger dataset. 



The spectra of IfD 98088 exhibit a complex pattern of two 
similar absorption line systems, with velocity separation 
changing from zero to ~200 kms"'^ and a period on the 
time scale of 5.9 days. In addition, one of the components 
has significant intrinsic line profile variability. In this situ- 
ation a spectral disentangling procedure is essential to sep- 
arate this effect from variable line blending due to orbital 
motion of the binary components and to obtain high-quality 
average spectra for abundance analysis. 

We empl oy the direct spectra l decomposition technique 
discussed bv iFolsom et al.l (|201G| ). Our algorithm operates 
in an iterative fashion, as follows. We begin with a set of 
approximate radial velocities for each component and ini- 
tial guesses of their spectra. Then the contribution of the 
less luminous component B is subtracted from the observed 
spectra, and all spectra are shifted to the rest frame of com- 
ponent A. The spectra are interpolated on a standard wave- 
length grid and co-added, yielding a new approximation of 
the primary spectrum. The same procedure is then used to 
update the approximation of the secondary spectrum. This 
sequence of operations is repeated up to convergence. Once 
converged spectra for the two components have been found 
the second part of the algorithm begins. In this part we use a 
least-squares minimisation routine to derive improved radial 
velocities (if desired) and correct the continuum normalisa- 
tion of the individual spectra using low-degree polynomials 
(up to this point the program does not distinguish between 
lines and continuum) . The first part of the program calculat- 
ing disentangled spectra and the second part of the program 
fitting radial velocities are alternated until the changes in 
parameters from one iteration to the next are below given 
thresholds. 

This spectral disentangling was applied to overlapping 
45 A segments of the spectra from 4500 A to 6590 A. While 
this wavelength range includes the Hq and II/3 Balmer lines, 
disentangling does not produce useful results for these lines. 
This is because the radial velocity variation between the 
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Table 1. Summary of MuSiCoS observations. The UT date, the Julian Date, the Stokes parameter, and the peak signal to noise ratio 
are given. The orbital phase of the observation was calculated using the ephemeris described in Table [2] Since the system appears to be 
synchronized, the rotational phase is assumed to be equal to the orbital phase. For each observation, four subexposures of 300s duration 
were acquired. 



Date 


HJD 


Phase Stokes Peak 


Velocity (kms ^) 




(2,450,000+) 


Param S/N 


Primary Secondary 



Jan 18, 1999 


1197.617 


0.358 


V 


190 


-8.0 




2.2 








Jan 18, 1999 


1197.643 


0.362 


Q 


240 


-9.2 


± 


2.3 








Jan 18, 1999 


1197.674 


0.368 


u 


230 


-6.3 


± 


2.2 








Jan 23, 1999 


1202.590 


0.200 


V 


140 


52.5 


± 


2.6 


-84.7 


± 


1.5 


Jan 23, 1999 


1202.617 


0.205 


Q 


210 


53.7 


± 


2.7 


-85.9 


± 


1.5 


Jan 24, 1999 


1203.555 


0.364 


V 


180 


-8.9 


± 


2.2 








Jan 24, 1999 


1203.581 


0.368 


Q 


200 


-9.9 


± 


2.3 








Jan 24, 1999 


1203.613 


0.374 


u 


200 


-7.7 


± 


2.2 








Mar 05, 2000 


1609.524 


0.112 


V 


190 


77.5 


± 


2.2 


-118.7 


± 


0.9 


Mar 07, 2000 


1611.519 


0.450 


Q 


230 


-33.7 


± 


2.2 


34.1 


± 


0.5 


Mar 07, 2000 


1611.554 


0.456 


u 


240 


-35.3 


± 


2.3 


34.5 


± 


0.5 


Dec 05, 2001 


2249.685 


0.520 


V 


100 


-53.2 


± 


2.4 








Dec 05, 2001 


2249.710 


0.525 


Q 


160 


-54.0 


± 


2.4 








Dec 05, 2001 


2249.742 


0.530 


u 


130 


-52.5 


± 


2.4 








Dec 07, 2001 


2251.693 


0.860 


Q 


270 


-28.6 


± 


2.3 


29.1 


± 


1.6 


Dec 07, 2001 


2251.725 


0.866 


u 


250 


-25.5 


± 


2.4 








Dec 07, 2001 


2251.751 


0.870 


V 


180 


-23.3 


± 


2.5 








Dec 08, 2001 


2252.690 


0.029 


Q 


230 


65.8 


± 


2.4 


-101.7 


± 


1.8 


Dec 08, 2001 


2252.722 


0.035 


u 


230 


67.3 


± 


2.2 


-103.8 


± 


2.5 


Dec 08, 2001 


2252.748 


0.039 


V 


160 


68.8 


± 


2.5 


-107.0 


± 


1.4 


Dec 09, 2001 


2253.687 


0.198 


Q 


250 


54.3 


± 


2.8 


-86.0 


± 


2.6 


Dec 09, 2001 


2253.715 


0.203 


u 


240 


52.9 


± 


2.5 


-85.0 


± 


1.8 


Dec 09, 2001 


2253.750 


0.209 


V 


230 


50.6 


± 


2.8 


-84.4 


± 


1.1 


Dec 12, 2001 


2256.710 


0.710 


Q 


270 


-71.5 


± 


2.6 


83.0 


± 


5.6 


Dec 12, 2001 


2256.741 


0.715 


u 


250 


-71.3 


± 


2.6 


82.5 


± 


2.8 


Dec 12, 2001 


2256.762 


0.719 


V 


110 


-71.2 


± 


2.6 


82.0 


± 


1.8 


Jan 05, 2002 


2280.694 


0.772 


Q 


230 


-63.2 


± 


2.5 








Jan 05, 2002 


2280.731 


0.778 


u 


200 


-61.3 


± 


2.4 








Jan 05, 2002 


2280.759 


0.783 


V 


130 


-60.5 


± 


2.4 









components of HD 98088 is insufficient to disentangle such 
broad lines, and because automatic continuum normalisa- 
tion of these broad features can be inaccurate. The radial 
velocities used were measured from the LSD profiles, pre- 
sented in in Table [T] and were not refined in the disentan- 
gling process due to intrinsic variability of the line profiles 
in the primary. 



Disentangling yields separate averaged high-quality 
{S/N ~ 1000) spectra of components A and B, as well as 
standard deviation curves which characterise the remaining 
discrepancy between the observations and composite model 
spectra. The standard deviation can be examined in the rest 
frame of the primary and secondary separately. A coher- 
ent excess in the standard deviation corresponding to a cer- 
tain absorption feature provides a straightforward, objective 
identification of intrinsic variability in that line. Performing 
this comparison, we found that all cases of intrinsic spectrum 
variability were associated with the lines of HD 98088 A. 
Variability in the primary was found in lines of Cr, Ba, Nd, 
and Eu, with weaker variability in lines of Mn and Sc. No 
clear variability was found in other lines, at the level of the 
noise in our spectra. 



5 FUNDAMENTAL PARAMETERS 

The luminosity ratio of HD 98088 was estimated using the 
Ha Balmer line cores. While the Balmer lines wings of the 
two components of HD 98088 are heavily blended, the cores 
are clearly separated in most of our observations. The depths 
of the Ha cores are only weakly dependent on tempera- 
ture, and thus can be used to estimate the fiux ratio of 
HD 98088 A to HD 98088 B at 6560 A. 

To model the H a line cores, we beg an with the T^b and 
log g estimates from lCarrier et all (|2002) . Rather than using 
synthetic Ha profiles, which can have inaccurate cores due to 
significant non-LTE effects, we used observed Ha profiles of 
the single stars HD 1086 51 and HD 73709. The stars were ob- 
served by Sho rlin et al.l (f2002) with MuSiCoS and have sim- 
ilar parameters to HD 98088 A (HD 108651: Tefi = 8100 K, 
log g = 4.36, vsini = 23 kms~^ ; HD 73709: TeS = 8080 K, 
log g = 4.02, sin i = 22 kms~^ IShorlin et al.ll2002l ). Obser- 
vations of the two stars were averaged to produce our model 
Ha line core. A composite binary Ha line was constructed by 
adding two model Ha lines, appropriately Doppler shifted, 
weighted by the flux ratio at 6560 A. This was then fit to 
the observed Ha line in the binary spectrum by varying the 
fiux ratio. This was repeated for seven observations, pro- 
ducing values of Fa/Fb of 3.0 to 3.3. An uncertainty was 
estimated by the range of values allowed by the noise in our 



© 2013 RAS, MNRAS OOO.nifTTl 



Orbit, composition, and magnetic field of HD 98088 5 




Figure 1. Measured radial velocities for the primary (full symbols) and secondary (empty symbols), phased with the orbital period of 
5.9051102 ± 0.0000023 days. Also shown is the best fit orbital curve to those velocities for the primary (solid line) and secondary (dashed 
line), using the ephemeris from Table [2] 




J 1 I 1 I 1 I 1 I 1 I 1 I i_ 

655.9 656 656.1 656.2 656.3 656.4 656.5 

Wavelength (nm) 

Figure 2. Sample fits to the observed combined Ha core of HD 
98088, for our best fit flux ratio at 6560 A (FR) of 3.1, and for 
two other flux ratios at ±2(t. 



observations, and our final best value is Fa/ Fb = 3.1 ± 0.2. 
An example of a fit to the Ha line core of an observation is 
given in Fig. (2] 

The flux ratio at 6560 A was then converted into a 
(bolometric) luminosity ratio by comparing with synthetic 
ATLAS9 flux distributions (|Kuruc j [l99^ ) for the atmo- 
spheric parameters in Table [3] at 6560 A, and assuming 
Lk/Lb = (reft,A/rcff,B)''(i?A/i?B)^ This was done itera- 



tively with the determination of effective temperature and 
surface gravity, discussed below. Using our flnal T^a, 8300 K 
for the primary and 7500 K for the secondary, this produces 
LA/iB = 3.62 ± 0.72 and R^/R^ = 1.55 ± 0.05. 

Effective temperatures for both components of 
HD 98088 were determined from the spectral energy dis- 
tribution (SED) of the system, since the Balmer lines were 
heavily blended. Since the system is not spatially resolved, 
simple photometry of the system is insufficient to determine 
the temperatures of both stars. 

Optical spectrophotometry was taken from the 
database of lAdelman et ah! (| 19891 ). Photometry in the UV 
was based o n observations f r om t he TDl satellite, as 
presented by [Thompson et al] (|l978l ) . IR photometry was 
take n from the Two Micron AU Sky Survey (2MASS) 
(Sk rutskie et al.ll2006l ). Synthetic model fluxes were taken 
from ATLAS9 models. Enhanced abundances (lOx solar) 
were used to reflect the chemically peculiar nature of the 
primary and, as we will show, the secondary. Both the mod- 
els and the observations were normalised to the fiux at 5000 
A, and are presented in magnitudes. The model fluxes were 
added, weighted by the luminosity ratio. 

Since HD 98088 is nearby (129.5 ± 6.7 pc) redden- 
ing of system is likely very small, and has t hus been ne - 
glected in our analysis. The reddening maps of lLuckel (| 19781 ) 
suggests that E(B — V) < 0.02. The detailed maps of 
ISchlegel et al.l l| 19981 ). which include all Galactic reddening 
along a line of sight, puts an upper limit on the reddening 
of E{B -V) < 0.05, though the true value for the HD 98088 
system is likely much less. Thus, following the recommenda- 
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tion of lLipski fc StepienI (|2008l ). since the reddening is small 
we neglect it. 

Fitting the SED provided well constrained effective tem- 
peratures for both stars. Since the primary dominates the 
flux in the UV, and the secondary's flux nearly matches 
the primary in the IR, precise T^s values for both stars 
can be derived. However, surface gravity cannot be deter- 
mined simultaneously for both stars with the same precision, 
since log g predominately affects the Balmer jump where 
flux from both stars is important. The luminosity weighted 
mean log g is 4.0, but there is a mild degeneracy between 
the two stars. The fit to the observed SED is presented in 
Fig. |3l and the best fit parameters for the primary are T^s 
= 8300 ±300 K and log g = 4.0 ±0.5, and for the secondary 
Teff = 7500 ± 300K and log g = 4.0 ± 0.5. 

HD 98088 was observed with Hipparcos, allowing us 
to determin e accurate luminosit ies. Based on the new re- 
duction by Ivan LeeuwenI |2003), the system has a par- 
allax of 7.72 ± 0.4 mas and a V magnitude of 6.173 ± 
0.007. The uncertainty in magnitude is the standard devia- 
tion of the Hipparcos measurements rather than the (sub- 
stantially smaller) quoted formal uncertainty, and agrees 
with the b and y band photometric variability found by 
ICatalano fc Leong (Il994|). The bol ometric correction for Ap 
stars from Landstreet et al.] l|2007h was used, since we find 
both components of HD 98088 have strong peculiarities, 
producing BCa = 0.08 and BCb ~ 0.09 (we adopted 
Af boi.a = 4.74). For comparison, the bolometric correction 
of iBalonal (| 19941 ) produces BCa = 0.01 and BCb = 0.04 . 
Following the recommendation of iLandstreet et"aL I (|2007l ). 
we adopted a conservative uncertainty on the bolometric 
correction of 0.1 in order to account for any potential sys- 
tematic uncertainties. Using the luminosity ratio we calcu- 
lated the bolometric magnitude difference between the pri- 
mary and the total system (rriA — mtot = 0.265 ±0.053), and 
the difference between the secondary and the total system 
(ttib— mtot = 1.662±0.191). The final luminosities we derive 
are shown in Table [3] 

With luminosities for both components, we can place 
the stars on the H- R diagram. By comp arison with the 
evolutionary tracks of ISchaller et al l (|l992l ) computed with 
Z = 0.02 and standard mass-loss, we determined masses for 
both components and an age for the system. The H-R di- 
agram positions of both stars are consistent with a single 
isochrone, suggesting that the stars are coeval. The age is 
determined from the H-R diagram position of the primary; 
since it is the more evolved star it allows us to determine a 
more precise age. The luminosities, masses, and age of the 
system are reported in Table [3] and the H-R diagram of the 
system is presented in Fig. |4l 

A few additional useful parameters can be calculated, 
and a few consistency checks can be made between our H- 
R diagram parameters and our dynamic parameters from 
the radial velocity curve. The mass ratio from the H-R di- 
agram Ma/Mb = 1.31 ± 0.08 agrees well with the more 
precise dynamic value Ma/Mb = 1.374 ± 0.008. Combining 
the dynamic Ma s\v? i with the H-R diagram Ma we derive 
the inclination of the orbital axis with respect to the line 
of sight iorb = 69 ± 2° . Using the secondary's mass we find 
the consistent value iorb = 67 ± 2°. Then the semi-major 
axis of the primary's orbit can be calculated, using the dy- 
namic OA sin i, yielding a a ~ 9.03 ± 0.10 Rq. Similarly, for 



Table 3. Summary of physical parameters for the HD 98088 sys- 
tem. 





HD 98088 A 


HD 98088 B 


d (pc) 


129.5 ±6.7 


iorb (°) 


69 


±2 


log age 


8.81 ±0.06 


Teff(K) 


8300 ± 300 


7500 ± 300 


log 9 (cgs) 


3.90 ±0.09 


4.16 ±0.12 


log L/Lq 


1.51 ±0.06 


0.95 ±0.10 


M (Mq) 


2.19 ±0.07 


1.67 ±0.08 


R (Rq) 


2.76 ±0.28 


1.77 ±0.24 



the secondary ae = 12 AO ± 0.15 Rq. We also calculate the 
inclination of the rotational axis of the primary, using the 
rotational period from the magnetic analysis and v sin i from 
the abundance analysis, and find irot,A = 68^^^°. While this 
is rather uncertain, it is consistent with the orbital inclina- 
tion, which one would expect for a tidally locked system. 
We can also calculate log g based on the star's H-R diagram 
masses and radii, finding log g = 3.90 ±0.09 for the primary 
and log g = 4.16 ± 0.12 for the secondary. These are consis- 
tent with the values from the SED, but significantly more 
precise. 

We can a l so com pare our values to the detailed work of 
[Carrier et all l|2002l '). and we find our physical parameters 
are fully consistent with theirs. Our orbital periods and in- 
clinations are consistent, as are our Tcff , log g, luminosities, 
masses and radii. 



6 CHEMICAL ABUNDANCES 

A detailed abundance analysis was performed for both 
HD 98088 A and B, using the disentangled spectra. Two in- 
dependent procedures were adopted, and independent analy- 
ses were performed by two of us. One analysis (performed by 
K.L.) employed th e LTE spectrum synthesis code Synth3 
(|Kochukhovfl2007l ) coupled to the IDL visualisation script 
BinMag. The other analysis (performed by C.P.F.) em- 
ployed a modified version o f the Zeeman spe ctrum synthe- 
sis code (|Landstreej|l988l : IWade et all 120011 ) which solves 
the polarized radiative transfer equations, again assuming 
LTE. Optimizations to the code for negligible magnetic fields 
were included, and a Levenberg-Marquardt minimisation 
routine was used t o aid in fitting the observed spectrum 
(|Folsom et al.ll2012h . 

In both cases, atomic data was ex tracted from the V i- 
enna Atomic Line Database (VALD) (jKupka et al.lll999l ). 
using an 'extract stellar' request, with temperatures tailored 
to the components of HD 98088. The model atmospheres 
used in the abun dance analysis were calculated using the 
LLMODELS code ijShulvak et al.ll2004l ). which produces LTE 
model atmospheres with detailed line blanketing. Prior to 
abundance analysis, the average disentangled spectra of the 
primary and secondary were corrected f or the continuum di - 
lution by the other component (see e.g. iFolsom et al.ll2008h 
using the theoretical, wavelength dependent light ratio pre- 
dicted by the adopted LLModels atmospheres. 

The Synth3 analysis involved directly fitting synthetic 
spectra to the deconvolved spectra of both components of 
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Figure 3. Fits to the spcctrophotomctric data using model binary flux distributions. The models are labelled by the T^ff of the primary 
and secondary, all for log g= 4.0 and a luminosity ratio of 3.6. The error bars in wavelength represent the full width of the filters for the 
UV and IR photometry. 




log(T.,) 

Figure 4. The positions of HD 98088 A and B on the H-R diagram. Evolutionary tracks (black) are labelled in unit s of solar mass e s, and 
isochrones correspond to the age of the system (8.81 it 0.06 log yr). Evolutionary tracks and isochrones are from ISchaller et al] lll992l 1 
for standard mass loss and Z = 0.02. 
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HD 98088. The fitting was performed independently on a 
number of small spectral windows. For the primary, 20 spec- 
tral windows were used with typical lengths of 100 A (rang- 
ing from 50 to 150 A long), running from 4500 to 6500 A, 
excluding Balmer lines. For the secondary, 13 windows were 
used typically with 150 A lengths, again running from 4500 
to 6500 A excluding Balmer lines. The fitting was performed 
by manually adjusting input chemical abundances, vsini, 
and microturbulence, then evaluating the fit of the synthetic 
spectrum to the observation with the BinMag visualisation 
tool. The best fit values from the individual windows were 
averaged to produce the final best values, and the uncertain- 
ties on the final best values were taken to be the standard de- 
viation of the values from individual windows. For elements 
with abundances derived only in three or fewer windows, the 
uncertainty was taken to include the line-to-line scatter as 
well as the range of abundances produced by varying T^s 
and log g by la. The final best fit values are reported in 
Table H 

The Zeeman analysis proceeded by fitting the spectra 
of both stars independently for chemical abundances, as well 
as vsini and microturbulence. Fitting was performed for 6 
independent spectral regions, between 200 and 400 A long, 
covering a spectral range from 4900 to 6500 A. The final 
best fit values adopted are averages of the results over all 6 
windows, with the adopted uncertainty as the standard devi- 
ation of those results. For elements with less than four lines 
available, the the uncertainty was estimated by eye instead, 
including potential normalisation errors, blended lines, and 
the scatter between lines. The final averaged best fit val- 
ues and uncertainties are presented in Table |4] The best fit 
ch emical abundan ces are plotted, relative to the solar values 
of lAsplund et al.l (,2009,) , in Fig. [S] for both components of 
HD 98088. 

6.1 Abundances of HD 98088 A 

The best fit synthetic spectra from both the Zeeman and 
the SynthS analyses generally reproduce the observed dis- 
entangled spectrum well. Samples of typical fits to the ob- 
servation for both analyses of HD 98088 A are presented 
in Fig. El Small differences between both synthetic spectra 
and the observation do exist. This is likely due to chemical 
stratification, which is common in Ap stars, or to blending 
with unknown rare-earth lines, which can also occur in Ap 
stars at these temperatures. 

HD 98088 A clearly shows typical Ap peculiarities, with 
strong overabundance of Si, Ca, Ti, Mn, Fe, and Ni, a dra- 
matic overabundance of Cr, and very strong overabundances 
of Ba, Ce, Pr, Nd, Sm and Eu. Ni and Sc have nearly solar 
abundances, as do C, O, Na, and Mg. 

In order to approximate the desaturation effect of the 
magnetic field of HD 98088 A in our spectra, a microtur- 
bulence of 2 kms^^ was assumed for the Zeeman analysis 
and a value of 1 kms~^ was used for the SynthS analy- 
sis. While the difference in assumed microturbulence was a 
consequence of the independent analyses, it serves to illus- 
trate the minor impact that small variations in microtur- 
bulence have on our results. There are a half-dozen careful 
analy ses of sharp lined Ap s tars published in the literature 
(e.g. iRvabchikova et a"l1ll997, : Kochukhoy et al. 200^, from 
Ryabchikova private communication), all of which find mi- 



croturbulence consistent with zero. Thus, while the typi- 
cal microturbulence in an Ap star is unknown, it is likely 
small. The use of microturbulence for HD 98088 A is simply 
to approximate line broadening due to the magnetic field, 
thereby reducing the computation time for spectra dramati- 
cally. Replacing the microturbulence with a dipole magnetic 
field like that determined in Sect. [7] would modify most of 
our abundances by 0.1-0.2 dex, which is generally within our 
uncertainties. 

The abundance of HD 98088 A derived from the two 
analyses agree very well. The abundances all agree within 
la, except for Ti and Mn at 1.2a and O at 2a. Both O 
abundances are based on one weak, heavily blended triplet 
of lines (6156-6158 A), thus an additional systematic effect 
may be present, or the uncertainties may be underestimated. 
The vsmi values agree at 1.4cr. No systematic differences 
between the results are found. 

An attempt was made to determine a spectroscopic T^a 
and log g by fitting metal lines of varying excitation poten- 
tials and ionisation stat es using the Zeema n program and 
the method described bv lFolsom et al.l l|2012l ). This attempt 
was unsuccessful. While we found log p = 4.1 ± 0.2 which 
appears to be accurate, we also found T^s = 9000 ± 200 
K which is clearly inconsistent with the observed SED. A 
likely cause of this discrepancy is chemical stratification in 
HD 98088 A, which could cause lines formed deeper in the 
stellar atmosphere to be anomalously strong, thereby pro- 
ducing an incorrect Teff estimate when metallic lines are 
used. 

The presence of chemical stratification in HD 98088 A 
was investigated using the Zeeman spectrum synthesis pro- 
gram. While stratification is not the focus of this study, these 
results could prove instructive for future studies, and they 
help explain the initial inaccurate spectroscopic T^s- Strati- 
fication was modelled with a simple step function, consisting 
of one abundance deeper in the atmosphere, another abun- 
dance higher in the atmosphere, and an optical depth for 
the transition. In this pilot study Fe stratification was in- 
vestigated, while other elements were assumed to have ho- 
mogeneous abundances. The optical depth scale used was 
the optical depth (due to continuous opacity) at the mid- 
dle of the wavelength region modelled. The Fe stratifica- 
tion parameters were added to the minimisation routine, 
and they as well as all other abundances and «sini were 
fit simultaneously, while the atmospheric parameters were 
held fixed at the values in Table [3] The fitting was per- 
formed for the 4900-5500 A, 5500-5700 A, and 6100-6500 
A regions independently. The resulting fits provide a clear 
improvement over a homogeneous model, both by visual 
inspection and by . Averaging over the three windows, 
and taking the range of values as the uncertainty, the Fe 
abundance higher in the atmosphere (lower optical depth) 
is Fe^pper = —4.5 ± 0.2 dex, the Fe abundance deeper in 
the atmosphere is Fcio^er = —2.2 ± 0.3 dex, and the tran- 
sition occurs at logr = —0.3 ± 0.2. While this represents 
strong Fe stratification, it is similar to what is seen in other 
Ap stars of similar Teg fe.g. IRvabchikova et al1l2002l . [2OO5I : 
iKochukhov et aLlbOOgl ). 

We then attempted to determine a spectroscopic Tcff 
simultaneously with Fe stratification (relying principally on 
the ionisation balance for Tcs). This was done by using the 
minimisation routine to fit Toff, Fe stratification, vsini, 
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Figure 5. Abundances relative to solar for HD 98088 A (top frame) and HD 98088 B (bottom frame), averaged over all spectral windows 
modelled. The results from bo t h the analysis with Zeeman (circles, black) and SynthS (squares, green) are presented. Solar abundances 
are taken from lAsplund et al.l ([200^. Arrows represent upper limits only. 



and other chemical abundances simultaneously, again using 
the 4900-5500 A, 5500-5700 A, and 6100-6500 A spectral re- 
gions. We find %ft = 8500±300 K, Fe^pper = -4.0±0.2 dex, 
Feio^er = -2.2 ± 0.4 dex, and logr = -0.5 ± 0.3. This T^s 
is consistent with the SED of HD 98088. The stratification 
result is is also consistent with the previous values, though 
the small difference in Fe„pper suggests it may be sensitive 
to the exact Tcfr used. As a note of caution, since only the 
stratification of Fe was modelled here, stratification of other 
elements could be influencing the derived Tcff potentially 
making it slightly inaccurate. However, this model does suc- 



ceed in reconciling the spectroscopic and spectrophotometric 
observations. 



6.2 Abundances of HD 98088 B 

We achieve good fits to the observed spectrum of 
HD 98088 B for both the SynthS and Zeeman analyses. 
Typical fits from both analyses of HD 98088 B are presented 
in Fig. [7] Both sets of best fit abundances for HD 98088 B are 
reported in Table |4j and plotted relative to solar in Fig. [5] 
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Table 4. Summary of abundances, Dsini, and microturbulence (5) for HD 98088 A and B, for both the analysis with SynthS and 
Zeeman. Chemical abundance are in units of log{Nx /Ntot)- For the Synth3 results, abundances fit manually with lines from three or 
fewer spectral regions are marked with an asterisk (*). F or the Zeeman results , abundances fit with three or less useful lines are marked 
with an double- asterisk (**). Solar abundances are from lAsplund et al. I l l2009t) . 





HD 98088 A 


HD 98088 B 


HD 98088 A 


HD 98088 B 


Solar 




SynthS 


SynthS 


Zeeman 


Zeeman 




V sin j (km s ^) 


22.0 ±1.5 


18.0 ±2.0 


24.4 ±0.8 


20.6 ± 0.8 




C (kms-i) 


1.0 ±0.5 


1.5 ± 0.5 


2.0 


1.9 ±0.2 





c 


-S 52 + 1 7* 


-4 62 + 40* 


-S 86+0 SO** 


—4 28 + 16 


-3.57 


o 


-4.21 ± 0.25* 


-S 61 + 20* 


-S S8+0 SO** 


— S 87 + 40** 

r.J*yJI _1_ U.'^U 


-3.31 


Na 


-fi 1 S + 54* 


-5 05 + 46 


-6 05+0 15** 


— 5 OS + SS** 


-5.76 


Mp- 


-4.14 di 0.23 


-5 02 + 50* 


-4 S2+0 06 


—4 78 + 26 

rt 1 J u _i_ yj . ^yj 


-4.40 




Q QO _|_ n 99 


<C-4 54 


Q Q9-|-n 9Q 


A AA -k n 1 9 


-4 49 


s 




^-4.7S 




-4.34 ± 0.24 


-4.88 


Ca 


-4.66 ± 0.19 


-6.68 ± 0.20 


-4.S9±0.S0 


-6.26 ±0.32 


-5.66 


Sc 


-9.10 ± 0.2S 


-9.70 ± 0.57 


-8.75±0.47 


-9.81 ±0.27 


-8.85 


Ti 


-6.56 ± 0.14 


-6.95 ± 0.15 


-6.24±0.22 


-7.05 ±0.13 


-7.05 


V 








-7.32 ±0.12 


-8.07 


Cr 


-4.05 ± 0.14 


-5.90 ± O.Sl 


-4.06±0.1S 


-5.48 ±0.18 


-6.36 


Mn 


-6.15 ± 0.25* 


-6.20 ± O.SO* 


-5.5S±0.46 


-5.92 ±0.19 


-6.57 


Fe 


-S.55 ± 0.2S 


-S.9S ± 0.11 


-S.59±0.04 


-3.95 ± 0.15 


-4.50 


Co 








-6.11 ±0.30** 


-7.01 


Ni 


-5.59 ± 0.35 


-5.19 ± 0.41 


-6.14±0.53 


-5.11 ± 0.18 


-5.78 


Cu 








-7.45 ±0.30** 


-7.81 


Zn 




-6.94 ± 0.50* 






-7.44 


Y 




-8.58 ± 0.72 


-9.46±0.S8 


-8.67 ±0.16** 


-9.79 


Zr 


-8.24 ± 0.55* 






-9.22 ± 0.30 


-9.42 


Ba 


-7.76 ± 0.10 


-7.46 ± 0.S7 


-7.8S±0.15 


-7.74 ±0.35** 


-9.82 


La 




-9.S7 ± 0.40* 




-9.26 ± 0.16 


-10.90 


Ce 






-9.0S±0.24 


-8.95 ± 0.15 


-10.42 


Pr 






-9.56±0.S9 


-9.63 ± 0.18 


-11.28 


Nd 


-9.40 ± 0.84 




-9.S4±0.20 


-9.12 ± 0.28 


-10.58 


Sm 






-8.85±0.17 


-9.13 ±0.14 


-11.04 


Eu 






-8.84±0.S1** 


-9.38 ±0.30** 


-11.48 



The abundances we find reveal that HD 98088 B is a 
clear Am star. We find strong overabundances of Cr, Mn, 
Fe, Co, Ni, and Y, and very strong overabundances of Ba, 
La, Cr, Pr, Nd, Sm and Eu. Sc and Ca are strongly under- 
abundant, C appears to be underabundant, and O and Mg 
are marginally underabundant. This pattern of abundances, 
with overabundant Fe-peak elements and rare-earths, and 
underabundant Ca and Sc is characteristic of Am stars. We 
also find a small microturbulence of « 2 kms^^. 

In the Zeeman analysis, microturbulence was deter- 
mined by including it as a free parameter in the analy- 
sis and fit simultaneously with chemical abundances. In the 
SynthS analysis, microturbulence was determined by fitting 
the synthetic spectrum to the observation by eye, together 
with chemical abundances. Microturbulenc e is often elevated 
in Am stars (typically 2 to 4 kms~^, e.g. iLandstreet et al.l 
l2009l ). and thus can have a significant impact on the abun- 
dance analysis of these stars. 

The results for the two abundance analyses agree very 
well. The chemical abundances, vsini, and microturbulence 
all agree within la, except for Ca and Cr which agree at 
l.lcr and 1.2a respectively. This confirms the accuracy of our 
methodology. Our uncertainties, if anything, may be slightly 



overestimated, since statistically one would expect a few of 
our results to deviate by more than la. 



We attempted to determine a spectroscopic T^b for 
HD 98088 B using fits to metallic lines, similarly to 
HD 98088 A. This was done using Zeeman and independent 
fits for all 6 windows used for the abundance analysis. The 
produced a value of T^g = 8000 ± 300 K, which is consistent 
with the result from the SED at the 1.2cr level. We prefer the 
SED value, since the spectroscopic T^g could be infiuenced 
by any small systematic errors in the disentangled spectrum, 
or by inaccurate or missing atomic data, particularly for the 
strongly overabundant rare-earth elements. 



This is the first report of chemical peculiarities in 
HD 98088 B, and a very rare case of a binary system with 
an Ap primary and an Am secondary. However, the ap- 
pearance of Am peculiarities in HD 98088 B may refiect 
the generally higher in cidence of Am peculiarities in low 
11 sin i stars dAbt fc Morrell 1995i ). and in binary systems 
(e.g. ICarauillat fc Prieuij|2o57l ). 
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Figure 6. Typical fits to the spectrum of HD 98088 A for both the Zeeman and SynthS results. Points represent the observation, the 
lighter grey (green) line is the Zeeman, and the darker grey (blue) line is the SynthS best fit spectrum. The dashed line is a synthetic 
spectrum with solar abundances. Elements which are major contributors to each line have been labelled. 
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Figure 7. Typical fits to the spectrum of HD 98088 B for both the Zeeman and SynthS results. Points represent the observation, the 
lighter grey (green) line is the Zeeman, and the darker grey (blue) line is the SynthS best fit spectrum. The dashed line is a synthetic 
spectrum with solar abundances. Elements which are major contributors to each line have been labelled. 
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7 MAGNETIC FIELDS OF THE 
COMPONENTS 

The magnetic fields in the photospheres of HD 98088 A and 
B were studied using the circular and linear polarisation 
within the LSD profiles, produced as a consequence of the 
longitudinal and transverse Zeeman effect. 

Fig. [8] illustrates the extracted Stokes /, V, Q and U 
LSD profiles, ordered according to orbital phase. In the left- 
most frame, the mean Stokes I profiles of both the primary 
and secondary stars are visible, and exhibit systematic dis- 
placement as a consequence of their orbital motion. The rel- 
ative weakness of the secondary's mean line is notable. In 
the remaining 3 panels, it is clear that Stokes V, Q and U 
Zeeman signatures are detected only at the position of the 
mean line of the primary star, indicating (at least at first 
glance) that the secondary is non-magnetic. Indeed, eval- 
uating the de tection probability according to the criteria 
introduced bv lDonati et al.l (Il997l ) indicates that significant 
Stokes V signal is detected in the mean line of the primary 
star at all phases, whereas no significant signal is ever de- 
tected in the mean line of the secondary. Finally, we point 
out the excellent phasing of the LSD profiles according to 
the orbital period. Although the data span more than 180 
orbits of the system, there is excellent fidelity of the detailed 
profile shapes and positions inferred to correspond to sim- 
ilar phases, and acquired both close together or far apart 
in time (e.g. phases 0.200/0.209 and 0.358/364). This is the 
most visually striking evidence for precise synchronization of 
the orbital period and the rotational period of the magnetic 
primary star. 

To quantify the magnetic properties of the two stars, we 
measured the longitudinal magnetic field Be of both compo- 
ne nts from the L S D pro files using the first-moment method 
of iRees fc Semell (|l979l ). We integrated the I/I^ and V/Ic 
profiles about their centr es-of-gravity vp in ve locity v, in the 
m anner implemented bv lDonati et al.l |l993) and corrected 
bv lWade et all l|2000l ): 

{v — Vo)V{v) dv 
Bi = -2.14 X 10" J . (1) 

Xzc [1- I{v)] dv 

In Eq. [U V{v) and I{v) are the V/L and I/L LSD pro- 
files, respectively. The wavelength A is expressed in nm and 
the longitudinal field Be is in gauss. The wavelength and 
Lande factor z correspond to the mean values employed in 
computing the LSD profiles, reported in Sect. (2] 

For the primary, the velocity integration limits for the 
evaluation of Eq.[T]were selected manually to include the en- 
tire observed Stokes V profile. In the case of the secondary, 
the limits were chosen to include the entire Stokes / profile. 
The centre-of-gravity for each component was also calcu- 
lated using these integration limits. Table [S] summarises the 
longitudinal magnetic field measurements for each date for 
which Stokes V observations were obtained. At phases at 
which the profiles of the primary and secondary star were 
blended, we measured only the longitudinal field of the pri- 
mary star by assuming that the secondary made no contribu- 
tion to the Stokes V profile, and adjusting the denominator 
of Eq. [1] by subtracting the average equivalent width of the 
secondary's LSD profile. Standard uncertainty propagation 



was used to determine the la error bars on the longitudinal 
magnetic field measurements. Uncertainties range from 15 
to 35 G for the primary and from 78 to 280 G for the sec- 
ondary component. A significant longitudinal magnetic field 
is detected in the primary component at all phases. On the 
other hand, no significant magnetic field is ever detected in 
the secondary component. 

We point out that the value of the longitudinal field of 
a component of an SB2 system inferred using Eq. [T]is unaf- 
fected by the continuum of the companion. This is because 
it is the continuum normalised Stokes / and V parameters 
{I/Ic and V/Ic, respectively) that appear in the denomi- 
nator and numerator of Eq. [1] Nevertheless, the inferred 
uncertainty, derived from photon-noise error bars propa- 
gated through the LSD procedure, is larger than that for 
the single-star case. 

From these basic magnetic measurements, we conclude 
that the primary component of HD 98088 hosts a strong 
magnetic field with a peak disc-averaged longitudinal com- 
ponent of about 1100 G, whereas the secondary shows no 
detectable magnetic field. 

7.1 Rotational period and magnetic geometry of 
the primary star 

The Lomb-Scargle periodogram of the longitudinal field 
measurements of the primary star displays a clear excess 
of power between 5.6 and 6.0 d, with many possible solu- 
tions in this range. The peak closest to the orbital period 
corresponds to a variability period of 5.90451 d, but with 
rather large uncertainties (~ ±0.0015 d). If we combine our 
new longitudinal field measurements with those reported by 
iBabcock (1958), this yields a period of 5.90510 ± 0.00003 d. 
This rotational period is in excellent agreement with the or- 
bital period (5.9051102 ± 0.0000023; Table [l]). Given this 
correspondence and the better precision of the orbital pe- 
riod, hereinafter we assume perfect synchronization of the 
system and adopt the orbital period as the rotational period 
of both the primary and secondary stars. 

The phased lon gitudinal field m easurements are shown 
along with those of lBabco"c^ l|l958h in Fig. (9] The signifi- 
cant improvement in precision of the MuSiCoS data relative 
to Babcock's measurements is evident. Although the his- 
torical measurements are noisy, it appears that the new Be 
curve may be shifted to slightly ('^ 100 — 150 G) higher val- 
ues than the previous measurements. If real, we attribute 
this shift to differences in the modelling and analysis meth- 
ods employed. Such a difference is not surprising, since the 
older measurements were based on (somewhat non-linear) 
photographic data, with line centres estimated by eye. Bab- 
cock's observations were made blueward of H/3, while our 
observations are predominately to the red of that line, thus 
the two datasets are affected by limb darkening somewhat 
differently, which could also partly explain the discrepancy. 
Our new measurements show clear and systematic depar- 
tures from a pure sinusoidal variation. Whereas first-order 
and second-order harmonic expansions fit phase variation of 
the new data rather poorly (reduced of 16.1 and 18.6, 
respectively), a third-order expansion achieves an excellent 
fit (reduced of 0.6). While these departures from a sinu- 
soidal variation may diagnose non-dipolar characteristics of 
the stellar magnetic field, it is more likely that they reflect 
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Figure 8. LSD Stokes IVQU profiles extracted from the observations of HD 98088. Tlie wealcer line of the secondary is visible in the 
leftmost frame. Polarisation signatures are observed at the position of the mean line of the primary, and follow the orbital motion of that 
line. No polarisation is observed at the position of the mean line of the secondary. Phases (xlOOO) are labelled in the leftmost frame. 



the non-uniform distributions of chemical elements whose 
lines were used to measure the field. 

To determine the magnetic geometry of the primary, 
we took advantage of the Stokes Q and U LSD profiles to 
measure the net li near polarization (NLP), integrated acros s 
the LSD profiles (Wado ct al.' 2000; Silv ester etlo] l2012h . 
These data are reported in Table [5] and are compared to 
the broadband linear polarization ( BBLP ) measurements of 
HD 98088 rep orted bvlLerov et al.l |l995l ). in Fig.[Tni As de- 
scribed by e.g. lSilvester et all l|2012l 'l. while NLP and BBLP 
result from the same physical mechanism, the methods of 
measurements are sufficiently different that their intercom- 
parison requires that one data set be arbitrarily shifted and 
scaled relative to the other. The NLP values reported in 
Table [S] and illustrated in Fig. [TO] include this scaling and 
shifting. The BBLP and NLP phase variations are generally 
compatible, although some differences in detail are certainly 
present. Th ese results ar e consi stent with the general resu lts 
reported by I Wade et al] (|2000h and lSilvester et al.l (|2012l ). 

In combination with the longitudinal field, the 
NLP/BBLP measurements provide an unambiguous esti- 
mate of the rotational axis inclination i, magnetic obliq- 
uity /3 and polar field strength Bd, under the assump- 
tion of a pure d ipolar surface magnetic field configura- 
tion (|Lerov et al .l [1995 1. A first-order harmonic fit to Bab- 
cock's data yields B^,max ~ 1050 ± 175 G and Bf.min = 



-1280±175 G, giving r = B£,niax/-B^,miii = -0.82±0.22. An 
analogous fit to the MuSiCoS measurements yields -Bf ,max = 
1185±65 G and B^.^in = -1000±65 G, for r = -0.84±0.09. 
The broadband linear polarization measurements of Leroy et 
al. give s = —0.065 ±0.13, where s is a (unitless) parameter 
describing the harmonic conte nt of the NLP/BB LP phase 
variation, defined according to iLerov et al.l (Il995 ). For the 
NLP variation, this value is -1-0.005 ± 0.20. The combined 
values of r and s yield unique values for the angles i and /3. 



Using this method Leroy ct al. ( 1996) reported i = 85± 
5° and /3 = 80 ± 5° for HD 98088 A. For the combination 
of BBLP + Babcock, we find i = 75 ± 15° and ^ = 73° ± 
3°. For BBLP -f MuSiCoS Bi, we find i = 77 ± 10° and 
13 = 73° ± 3°. Finally, for NLP + MuSiCoS Be, we find 
i = 77± 10° and /3 = 76° ± 5° . All of these values are in good 
mutual agreement, indicating that the systematic differences 
between the various data sets are not a major source of error. 
We also note that the derived values of the inclination i are 
formally consistent with the orbital inclination determined 
from the masses (68 ±3°). 

If we adopt the orbital inclination as the inclination of 
the rotation axis of the primary, along with /3 = 75 ± 5° , we 
derive the surface polar strength of the primary's magnetic 
dipole to be 3850±450 G using the MuSiCoS measurements. 
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Figure 9. Phased longitudinal magnetic field measurements of 
HD 98088 A. Filled squares - MuSiCoS measurements. Open tri- 
angles - Measurements reported by Babcock (1958). The vertical 
dashed lines indicate the approximate quadrature phases of the 
radial velocities. 
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Figure 10. Net linear polarisation of HD 98088 A, compared to 
the br oadband linear polari zation measurements of iLerov et al.l 
(|l995l) . As discussed by e.g. ISilvester et al.l ||2012|'| . the MuSiCoS 
measurements have been arbitrarily shifted and scaled to repro- 
duce as closely as possible the broadband variations. 



7.2 Upper limit on magnetic field of the 
secondary star 

No magnetic signature is detected in tlie mean line of tiie 
secondary, nor is tiie longitudinal magnetic field detected 
witii significance. Ifowever, we can estimate an upper limit 
on the dipole component of the secondary. Assuming the 
secondary's rotation is synchronized with the orbit, i.e. with 
rotational period equal to Porb and rotation axis inclination 
i = 68°, the median longitudinal field error bar of 160 G 
translates into a So" upper limit on the surface dipole of 
roughly 1550 G. 



8 DISCUSSION AND CONCLUSIONS 



We find strong chemical peculiarities in HD 98088 A that 
are characteristic of a cool Ap star. In HD 98088 B we also 
find strong chemical peculiarities indicative of an Am star. 
We also find evidence of strong Fe stratification in the at- 
mosphere of the primary. The introduction of stratification 
was necessary to reconcile the spectroscopic and spectropho- 
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Figure 12. Schematic of the orbit of HD 98088 A (circle with 
arrow) from the reference frame of HD 98088 B (filled circle). A 
dipole magnetic field for HD 98088 A is indicated with an arrow, 
assuming the dipole is aligned with the secondary at periastron. 
The dotted line represents the line towards the observer. The 
stellar radii and orbital separation are approximately to scale. 
Phases where the observed longitudinal magnetic field (B;) is a 
maximum, minimum, and zero are shown, as are phases where 
the stars are aligned along the line of sight. Note that the phases 
where the radial velocities of the components are equal are shifted 
slightly from the phase where the stars are aligned along the line 
of sight (by +0.01 and -0.03 in phase), due to the eccentricity of 
the orbit. 



tometric temperatures of HD 98088 A. While Am stars are 
often found in close binary systems, HD 98088 represents one 
of very few close binary systems with an Ap component. 

Both components of HD 98088 have similar masses. Toff, 
log g, and vsmi, yet the components display very different 
chemical peculiarities. Fig. [11] illustrates the difference in 
the derived abundances between the primary and the sec- 
ondary. The components of HD 98088 are likely coeval, but 
a qualitative difference is the presence of a strong magnetic 
field in the primary. This suggests that the observed differ- 
ences in chemical abundances are due to the impact of the 
magnetic fiel d on atomic diffusion in the atmosphere of the 
primary (e.g. iMichaud et al1ll98ll '). 

We find a strong, predominately dipolar, magnetic field 
in HD 98088 A, with a longitudinal strength varying from 
-f 1170 to -920 G. Conversely, we find no magnetic field in 
the secondary with longitudinal field error bars down to 112 
G, and a typical error bar of 160 G. We also find no signal 
from the secondary in linear polarisation, while the magnetic 
field of the primary is detected in linear polarisation at all 
phases. The predominately dipolar nature of the magnetic 
field in the primary is shown by the longitudinal field and 
net linear polarisation curves, and by the simple structure 
of the V , Q, and U LSD profiles. 

It is not clear why the primary has a strong magnetic 
field and the secondary has no magnetic field when the two 
stars have similar masses and temperatures, and are appar- 
ently coeval with identical ages. Indeed, the origin of the 
magnetic fields in Ap stars is one of the enduring mysteries 
that these stars present. 
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Table 5. Magnetic measurements obtained from MuSiCoS spectra. All phases are computed from the orbital ephemeris. HJD is from 
the Stokes V observation. Typical range of phase corresponding to a single Stokes VQU sequence is approximately 0.01 X Po^b- 



Date 


HJD 


Orbital 






NLP 




(2,450,000+) 


Phase 


Primary 


Secondary 


Q/I (%) 


U/I (%) 


18 Jan 99 


1197.617 


0.358 


1048 ± 15 




0.015 ±0.004 


0.048 ± 0.004 


23 Jan 99 


1202.590 


0.200 


384 ± 33 


51 ± 166 


0.027 ±0.007 




24 Jan 99 


1203.555 


0.364 


1062 ± 19 




0.007 ±0.005 


0.043 ± 0.005 


05 Mar 00 


1609.524 


0.112 


-370 ± 24 


-210 ± 135 






07 Mar 00 


1611.519 


0.450 






0.004 ± 0.005 


0.067 ± 0.005 


05 Dec 01 


2249.685 


0.520 


1167 ±27 


-195 ± 193 


0.008 ± 0.008 


0.086 ± 0.009 


07 Dec 01 


2251.751 


0.870 


-916 ± 19 


53 ±78 


-0.035 ±0.004 


0.039 ± 0.005 


08 Dec 01 


2252.748 


0.039 


-762 ± 28 


-92 ± 133 


0.041 ± 0.006 


0.062 ± 0.006 


09 Dec 01 


2253.750 


0.209 


435 ± 21 


226 ± 112 


0.023 ± 0.006 


0.102 ±0.006 


12 Dec 01 


2256.762 


0.719 


-291 ±35 


423 ± 280 


-0.034 ±0.005 


0.094 ± 0.006 


05 Jan 02 


2280.759 


0.783 


-674 ± 30 


-159 ± 156 


-0.036 ±0.006 


0.070 ± 0.007 
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Figure 11. Differential comparison of the chemical abundances found for the primary and secondary of HD 98088. Values from both 
the Synth3 (squares, green) analysis and the Zeeman (circles, black) analysis are presented. 



We derived improved orbital parameters for the 
HD 98088 system, finding a period of 5.9051102 ±0.0000023 
days. We also derive a rotational period for HD 98088 A 
of 5.90510 ± 0.00003 days, which is consistent with the or- 
bital period. This provides strong evidence that the system 
is tidally locked, with the rotational and orbital periods syn- 
chronised. We derive an inclination of the orbital axis of the 
system of iorb = 69 ± 2°, which is consistent with the incli- 
nation of the rotational axis of the primary based on its ro- 
tational period, vsmi, and radius, as well as the inclination 
derived from the linear polarisation under the assumption 
of a dipole magnetic field. However, the system has a mod- 
est but significant eccentricity of 0.1840 ± 0.0025, which is 
clearly visible in the radial velocity curve presented in Fig.[T] 
Thus while the system is synchronised, the orbits apparently 
have not had the time to circularise. 

iHud (|l98ll ) showed that, in binary systems where orbital 
angular momentum dominates over rotational angular mo- 
mentum, tidal interactions cause the orbital and rotational 
axes to become parallel, and the orbital and rotational peri- 
ods to become synchronised, much more quickly than they 
circularise the orbits of the system. Using Eq. 22 of iHutI 



('l98l') (and a ssuming an approximate radius of gyration of 
0.2, following iStepiehll2000 ). the ratio of orbital to rotational 
angular momentum is ~150 for the pri mary, and ^500 for 
the secondary. This suggests (|Hutlll98ll . Eq. 33) that circu- 
larisation of the orbits would be slower by a factor of '^lOO 
than the synchronisation of orbital and rotational periods. 
Thus it is certainly possible for HD 98088 to be rotationally 
synchronised but still have elliptical orbits. However, we find 
an age for the system of 8.81 ± 0.06 logyr, which suggests 
that tidal interactions would have to be very inefficient for 
the system to still have some ellipticity. 

Modelling the magnetic field of HD 98088 A as a dipole, 
we find a dipole strength of 3850±450 G and an obliquity an- 
gle oi j3 = 75±5°. The longitudinal magnetic field reaches an 
extremum at very near the quadrature orbital phases (where 
the radial velocities of both components are near zero, see 
Figs. [T]and[9l. Thus the magnetic dipole is close to being 
aligned along the axis between the two stars, with the pos- 
itive magnetic pole always pointing near the secondary, as 
illustrated in Fig. 1121 This analysis does not account for the 
influence of abundance spots on the Bp curve, and if this 
effect were accounted for (with a self consistent modelling 
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procedure, such as magnetic Doppler imaging) the alignment 
could be exact. This is likely significant, though the inter- 
pretation of this observation is not immediately obvious. It 
could be that the tidal interaction between the components 
somehow influenced the stable magnetic field geometry of 
the primary. For example if the ma gnetic field had a foss il 
origin as in the models of .Braithwaite fc Nordlundl (|2006l ) , 
then tidal effects could affect the evolution of the magnetic 
field towards its stable configuration. However, magnetic 
fields in Ap stars appear t o originate very early in a star's 
lifetime (|Wade et alfbOOSl . in a few Myr), and the orbital 
parameters may have changed since then. One could specu- 
late that tidal interactions over the intervening years could 
have caused the field to evolve into its present configuration. 
Alternately, if the magnetic field perturbed the structure of 
the primary sufficiently to produce a small over-density at 
the magnetic poles, the system could settle into its current 
configuration through tidal effects. 

HD 98088 represents a very rare case of a magnetic Ap 
star in a tidally locked binary system with an Am star. The 
magnetic field in the Ap star is predominately dipolar, with 
the positive pole always pointing towards the secondary. 
Further studies could help clarify the interpretation of this 
result. Magnetic Doppler imaging of HD 98088 A would be 
particularly valuable, as this would provide a more detailed 
model of the magnetic field, as well as a map of chemical 
inhomogeneities across the surface of the star. These would 
be particularly interesting in the context of tidal interac- 
tions with the secondary. HD 98088 is an ideal target for 
the Binarity and Magnetic Interactions in Stars (BinaMIcS) 
project, and further observations and analysis are planned 
by the collaboration. 
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